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Abstract 

We demonstrate how to construct three-dimensional compact hyperbolic polyhedra using Newton's Method. 
Under the restriction that the dihedral angles are non-obtuse, Andreev's Theorem provides as necessary and 
sufficient conditions five classes of linear inequalities for the dihedral angles of a compact hyperbolic polyhedron 
realizing a given combinatorial structure C. Andreev's Theorem also shows that the resulting polyhedron is 
unique, up to hyperbolic isometry. Our construction uses Newton's method and a homotopy to explicitly follow 
the existence proof presented by Andreev, providing both a very clear illustration of proof of Andreev's Theorem 
as well as a convenient way to construct three-dimensional compact hyperbolic polyhedra having non-obtuse 
dihedral angles. 

As an application, we construct compact hyperbolic polyhedra having dihedral angles that are (proper) integer 
sub-multiples of vr, so that the group F generated by reflections in the faces is a discrete group of isometrics of 
hyperbolic space. The quotient H^/F is hence a compact hyperbolic 3-orbifold, of which we study the hyperbolic 
volume and spectrum of closed geodesic lengths using SnapPea I58| . One consequence is a volume estimate for a 
"hyperelliptic" manifold considered in I35| . 

1 Introduction 

Andreev's Theorem [21 El provides a complete characterization of compact hyperbolic polyhedra having 
non-obtuse dihedral angles. See also [^SJ^J for an alternative exposition on the classical proof. Other 
approaches to Andreev's Theorem can be found by Rivin and Hodgson |37[ 1211 . Thurston Marden 
and Rodin [30], and Bowers and Stephenson In this paper we show that the classical proof from 

[HI El ESI El is constructive when combined with Newton's Method for solving nonlinear equations. 

Combinatorial descriptions of hyperbolic polyhedra that are relevant to Andreev's Theorem fall into 
three classes, simple, truncated, and compound, all defined later in this section. The proof in [251 141) 
provides an explicit continuous path in the space of polyhedra deforming a given simple polyhedron P 
to one of two which are easily constructed by hand: the A^-faccd prism Ptn and the A^-faced split prism 
Dm- We use Newton's method to follow such a path backwards deforming a computer realization of 
Ptn or _Djv to a computer realization of the desired polyhedron P. This technique, which has been well 
studied in the literature, is known as the homotopy method [51 II II RTl I46L H51 149L 17] . We illustrate the 
construction of simple polyhedra in Sections 12.512.61 and 12.71 

A similar deformation, again using Newton's method, allows us to construct truncated polyhedra 
from simple polyhedra. We demonstrate this deformation in Section ITHl In Section [T^ wc show how to 
construct a compound polyhedron as a gluing of two appropriate truncated polyhedra. 

In this way, our program graphically illustrates Andreev's proof of existence for explicit examples. In 
fact, writing this program and working through Andreev's proof for some specific examples led to the 
detection an error in the proof of existence, which has been corrected in [251 141) . 
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A further benefit of this program is the construction of polyhedra whose dihedral angles are proper 
integer sub- multiples of tt. As a consequence of Poincares polyhedron theorem |20j . the group T generated 
by reflections in the faces of such a polyhedron is a discrete group of isometries of hyperbolic space. The 
quotient H'^/r is hence a compact hyperbolic 3-orbifold in which we study the hyperbolic volume and 
spectrum of closed geodesic lengths using SnapPea • Such orbifolds and their covering manifolds have 
been studied extensively [^Hl OHl ESI ^1 12H1 OH]' In fact the first example of a closed hyperbolic 
3- manifold was obtained in this way in 1931 by Lobell One consequence of our study is a volume 
estimate for a "hyperelliptic" manifold considered in |35) . 

The reader should note that there are already excellent computer programs for experimentation with 
hyperbolic 3-manifolds. The program SnapPea |HH1 constructs hyperbolic structures on knot and link 
compliments, as well as the hyperbolic Dehn surgeries on these compliments. SnapPea provides for the 
computation of a variety of geometry invariants of the computed hyperbolic structure. (See also 0].) The 
program Snap |18l I14| provides a way of computing arithmetic invariants of hyperbolic manifolds. Both 
of these programs are quite easy to use and have allowed for vast levels of experimentation, including a 
nice census of low-volume hyperbolic manifolds and orbifolds. 

The experimentation done in this paper with the hyperbolic orbifolds obtained from polyhedral re- 
flection groups is very modest in comparison. However, it is an alternative way to construct hyperbolic 
structures on certain orbifolds (and, in the future, possibly on manifold covers of these orbifolds) in a 
way that these structures can nicely be studied by SnapPea (as well as Snap, and other software, in the 
future). 

Let E^'^ be R"* with the indefinite metric ||x|p = —Xq + xl + X2 + x\. In this paper, we work in the 
hyperbolic space given by the component of the subset of E^'"^ given by 

I ||2 2 I 2 I 2 I 2 1 

|x|| = -Xq + a^i + 2:2 + 2^3 = -1 

having > 0, with the Riemannian metric induced by the indefinite metric 

^dxQ + dxl + dx\ + dx^. 

The hyper-plane orthogonal to a vector v G i?'^'^ intersects H"^ if and only if (v, v) > 0. Let v G E'^''^ 
be a vector with (v, v) > 0, and define 

Pv = {w e tf|(w^,v) = 0} and = {w e H'''|(w,v) < 0} 

to be the hyperbolic plane orthogonal to v and the corresponding closed half space. 

If one normalizes (v, v) = 1 and (w, w) = 1 the planes and P-^ in H"^ intersect in a line if and only 
if (v,w)^ < 1, in which case their dihedral angle is arccos(— (v, w)). They intersect in a single point at 
infinity if and only if (v,w)^ = 1; in this case their dihedral angle is 0. 

A hyperbolic polyhedron is an intersection 

n 

having non-empty interior. 

We will often use the Poincare ball model of hyperbolic space, given by the unit ball in with the 
metric 

dxf + dx\ -\- dxl 
(l-||xP)2 
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and the upper half-space model of hyperbolic space, given by the subset of R'' with > equipped with 
the metric 

dxj + dxl + dxl 



Both of these models are isomorphic to H"^. 

Hyperbolic planes in these models correspond to Euclidean hemispheres and Euclidean planes that 
intersect the boundary perpendicularly. Furthermore, these models are conformally correct, that is, the 
hyperbolic angle between a pair of such intersecting hyperbolic planes is exactly the Euclidean angle 
between the corresponding spheres or planes. 

1.1 Combinatorial polyhedra and Andreev's Theorem 

A compact hyperbolic polyhedron P is topologically a 3-dimensional ball, and its boundary a 2-sphere 
§^ . The face structure of P gives §^ the structure of a cell complex C whose faces correspond to the faces 
of P. 

Considering only hyperbolic polyhedra with non-obtuse dihedral angles simplifies the combinatorics 
of any such C: 

Proposition 1.1 (a) A vertex of a non-obtuse hyperbolic polyhedron P is the intersection of exactly 3 
faces. 

(b) For such a P, we can compute the angles of the faces in terms of the dihedral angles; these angles are 
also < 7r/2. 

See EH EH. 

The fundamental axioms of incidence place the following, obvious, further restrictions on the complex 

C: 

• Every edge of C belongs to exactly two faces. 

• A non-empty intersection of two faces is either an edge or a vertex. 

• Every face contains not fewer than three edges. 

Any trivalent cell complex C on that satisfies the three conditions above is an abstract polyhedron. 
Since C must be a trivalent cell complex on S^, its dual, C*, has only triangular faces and the three above 
conditions ensure that it is a simplicial complex on S^. The figure below shows an abstract polyhedron 
C drawn in the plane (i.e. with one of the faces corresponding to the region outside of the figure.) The 
dual complex is also shown, in dashed lines. 
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We call a simple closed curve T formed of k edges of C* a k- circuit and if all of the endpoints of the 
edges of C intersected by T are distinct, we call such a circuit a prismatic k- circuit. The figure below 
shows the same abstract polyhedron as above, except this time the prismatic 3-circuits are dashed, the 
prismatic 4-circuits are dotted, and the dual complex is not shown. 




We say that a combinatorial polyhedron C is simple if it has no prismatic 3-circuits, truncated if C 
has prismatic 3-circuits and each surrounds on one side a single triangular face, and otherwise we call C 
compound. The combinatorial polyhedron shown in the two above diagrams is compound. 

Theorem 1.2 Andreev's Theorem 

Let C be an abstract polyhedron with more than 4 faces and suppose that non-obtuse angles are 
given corresponding to edge Ci of C . There is a compact hyperbolic polyhedron P whose faces realize C 
with dihedral angle an at each edge e,; if and only if the following five conditions all hold: 

1. For each edge Ci, a.i > 0. 

2. Whenever 3 distinct edges e^, e^, eu meet at a vertex, a; + a^ + a^ > tt. 

3. Whenever V is a prismatic 3- circuit intersecting edges e^, Cj, Ck, a-; + a^ + a^ < tt. 

4. Whenever T is a prismatic ^-circuit intersecting edges e^, Cj, Ck, ei, then a.i + a.j + a^ + a; < 27r. 

5. Whenever there is a four sided face bounded by edges ei, 62, 63, and 64, enumerated successively, 
with edges 612,623,634,641 entering the four vertices (edge e^ connecting to the ends of Ci and Cj), 
then: 

ai + a3 + ai2 + a23 + a34 + a4i < 37r, and 
a2 + a4 + ai2 + a23 + a34 + a4i < 37r. 

Furthermore, this polyhedron is unique up to isometrics of . 

Corollary 1.3 If C is simple, i.e. has no prismatic 3-circuits, there exists a unique hyperbolic polyhedron 
realizing C with dihedral angles 2tt/5. 

For a given C, let E be the number of edges of C. The subset of (0, 7r/2]^ satisfying these linear 
inequalities will be called the Andreev Polytope, Ac. Since Ac is determined by linear inequalities, it is 
convex. 

Andreev's restriction to non-obtuse dihedral angles is emphatically necessary to ensure that Ac be 
convex. Without this restriction, the corresponding space of dihedral angles, Ac, of compact (or finite 
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volume) hyperbolic polyhedra realizing a given C is not convex |^. In fact, the recent work by Diaz 
[TB| provides a detailed analysis of this space of dihedral angles Ac for the class of abstract polyhedra C 
obtained from the tetrahedron by successively truncating vertices. Her work nicely illustrate the types of 
non-linear conditions that are necessary in a complete analysis of the larger space of dihedral angles Ac. 

The work of Rivin |39[ shows that the space of dihedral angles for ideal polyhedra forms a convex 
polytope, without the restriction to non-obtuse angles. (See also |19|.) 

Notice also that the hypothesis that the number of faces is greater than four is also necessary because 
the space of non-obtuse dihedral angles for compact tetrahedra is not convex ^^l- Conditions (1-5) 
remain necessary conditions for compact tetrahedra, but they are no longer sufficient. 

Bao and Bonahon ^0] prove a similar classification theorem for hyperideal polyhedra. Finally, the 
papers of Vinberg on discrete groups of refiections in hyperbolic space 1^ 1561 IS7] are also closely 
related, as well as the work of Bennett and Luo ^3 and Schlenker |43l I44L 1^ . 

Much attention has been focused on Andreev's Theorem from the viewpoint of circle packings and 
circle patterns. Given a polyhedron P in the upper half-space model of H'^, the planes supporting the 
faces of P intersect the boundary at infinity = in a pattern of circles (and straight lines) each with 
an orientation specifying "on which side" is the polyhedron P. Similarly, from such a pattern of circles 
and orientations one can re-construct a polyhedron P. 

The works of Thurston [SI] , Marden and Rodin PHI , and Bowers and Stephenson ^tU follow this 
approach to Andreev's Theorem. In fact, there is a beautiful computer program known as Circlepack |5U| . 
written by Ken Stephenson, that computes circle packings and patterns of circles with specified angles 
of overlap. All of the proofs from this point of view use the conformal structure of the Riemann sphere 
C = doo^^ and use the correspondence between conformal automorphisms of C with isometrics of H'^. 

Instead of using the conformal structure on SooEI^, in this paper we will work specifically with the 
metric structure of H^. (However, there is certainly some significant overlap with the results in |51II3()|[T^ 
and with the capabilities of the computer program CirclePack |5()).) 

We will now explain the implementation of a computer program whose input is the combinatorial 
polyhedron C and a dihedral angle vector a e Ac and whose output is a hyperbolic polyhedron realizing 
the pair (C, a). 

1.2 An example 

The following figure shows an explicit example of the data (C, a) and the resulting polyhedron displayed 
in the conformal ball model using the computer program Geomview 

1.3 Outline of the proof of Andreev's Theorem 

In this section, we recall the major steps from the proof of Andreev's Theorem that were presented in 

[23111]. 

Let C be a trivalent abstract polyhedron with N faces. We say that a hyperbolic polyhedron P C 
realizes C if there is a cellular homeomorphism from C to dP (i.e., a homeomorphism mapping faces of 
C to faces of P, edges of C to edges of P, and vertices of C to vertices of P). Wc will call each isotopy 
class of cellular homeomorphisms : C — > dP a marking on P. 

We defined Vc to be the set of pairs (P, (f>) so that (p is a. marking with the equivalence relation that 
(P, <j>) ~ (P', (f)') if there exists an automorphism p : H'^ ^ such that p(P) = P', and both <j)' and po 
represent the same marking on P'. 

Proposition 1.4 The space Vc is a manifold of dimension 3N — 6 (perhaps empty). 

The proof is relatively standard and can be found in |25l E] . 

Since the edge graph of C is trivalent, the number E of edges of C is the same as the dimension of 
Vc- Given any P g Vc, let a(P) = (ai, a2, as, ...) be the i5-tuple consisting of the dihedral angles of P 
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Andreev's Theorem 




Figure 1: 
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at each edge (according to some fixed numbering of the edges of C) . This map a is obviously continuous 
with respect to the topology on Vc^ which it inherits from its manifold structure. 

We let be the subset of Vc consisting of polyhedra with non-obtuse dihedral angles. To establish 
Andreev's Theorem, we proved the following statement: 

Theorem 1.5 For every abstract polyhedron C having more than four faces, the mapping a : — > Ac 
is a homeomorphism. 

There were two major steps: 

Proposition 1.6 IfV^ ^ 0, then a : ^ Ac is a homeomorphism. 

We checked that a{Vc) C Ac by showing that conditions (l)-(5) are necessary. There is an open 
subset Vc C Vc containing Vc on which one can prove that a : Vc ^ is injective, using a modification 
of Cauchy's rigidity for Euclidean polyhedra. This gives the uniqueness part of Andreev's Theorem. Using 
invariance of domain, it also gives that a : Vc — is a local homeomorphism. Because Vc C Vc, a 
restricted to Vc is a local homeomorphism, as well. 

We then showed that a : Vc — > Ac is proper, which amounts to showing that if a sequence of 
polyhedra P,; in Vc degenerate (i.e. leave Vc) then the sequence a{Pi) tends to dAc- The fact that 
a : Vc Ac is a proper local homeomorphism was sufficient to show that a{Vc) is open and closed in 
Ac. 

Proposition 1.7 If Ac ^ 0; then ^ 0. 

The second step was much more difficuh because for each C with non-empty Ac one needed to 
construct some polyhedron reahzing C (with non-obtuse dihedral angles). In fact the proof of Proposition 
Il.Tl outlines a scheme for how to construct a polyhedron realizing C . The remainder of this paper outlines 
how to follow this scheme explicitly on the computer using Newton's Method and a homotopy method. 



2 A method for constructing Andreev polyhedra 
2.1 Representing polyhedra on the computer 

All the constructions of polyhedra done in this paper are using Matlab ^ or the Free Software Foundation 
alternative Octave [21, and all of the polyhedra are displayed in Geomview When doing calculations, 
we represent a hyperbolic polyhedron P having faces by specifying A'^ outward pointing normal vectors 
vectors Vi, • • -Vtv each with (v^, v^) — 1, so that P = HiLi ^vi- 

Although such a list of N vectors is sufficient to specify P, in order to avoid repeated computation of 
the combinatorial structure of P from these vectors we additionally specify the adjacency matrix and a 
list of all plane triples meeting at a vertex. These three items are described in a Matlab struct P, with 
P.faces, Padjacency, and P.vert holding the data mentioned above, respectively. 

For example, the data for the polyhedron shown in Section fL^ is stored in Matlab as: 



[28x3 double] 
[16x4 double] 
[16x16 double] 



NeH_poly . f a( 
36.5078 

4.9237 
-0.0000 

4.5134 

2.7290 
13.3691 
65.0863 
35.9576 
51.5713 

5.8378 
-0.0000 



-10.7624 
-1.5291 
0.8660 
-2 . 1988 
-1.9854 
-5.0338 

-19.6363 
-9.6209 

-13.3145 
-0.5352 
0.0000 



-0.3090 
-1.2342 
-0.5000 
-2 . 3943 
-0.3091 
-0.3090 
-2.7939 
-1.5515 
-0.3090 
-0.3090 
1.0000 



-34.8983 
-4 . 6240 
-0.0000 
-3.2868 
-2.1000 
-12.4216 
-61.9987 
-34.6262 
-49.8320 
-5.8905 
0.0000 



1 . 2179 -1 . 4082 -0 . 7071 . 0000 

-7 . 8692 2 . 1329 3 . 6943 6 . 6879 

3 .4039 -1 . 5744 -2 . 7269 -1 . 6344 

-1 . 0781 . 5773 -0 . 0000 -1 . 3524 

-2 . 1544 . 9964 1 . 7260 -1 . 2921 
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3 4 14 
6 5 11 

4 12 5 
12 5 11 

16 7 

17 8 

18 9 
1 9 11 

10 9 11 

1 6 11 

2 3 13 



2 10 3 

11 10 3 

2 4 13 

3 4 13 

3 12 14 

4 12 14 
11 3 15 

11 12 15 
3 12 16 
3 15 16 

12 15 16 



We display the polyhedra in Geomvicw using the hyperbohc mode and specifying the conformal ball 
model. The file format most convenient for our use is the Object File Format, "New.poly.off." The first 
line of an Object format file specifies the number of vertices, the number of faces, and the number of 
edges of P in that order: num.vert num_faces num.edges. 

The next block of data is a list of the coordinates of vertices as points in the unit ball. (In fact, 
these are the coordinates of points in the Projective Model for H"^, not the Poincare ball model which we 
describe in the introduction.) The last block of data is a list of the faces with each face given by vertexi 
vertex2 ... vertex^ colorspec, where the faces is spanned by vertexi vertex2 ... vertex^ and colorspec is 
an integer telling Geomvicw what color to assign to the face. 



28 16 42 

0.093414 0.626297 -0.759378 
0.668701 -0.423986 0.508660 
0.480895 0.480927 -0.729094 
0.533074 -0.046431 -0.835831 
0.000602 -0.321164 0.944739 
-0.109909 -0.298793 0.946284 
-0.257482 -0.413626 0.871178 
-0.241511 -0.517345 0.817366 
-0.394737 -0.502851 0.762029 
0.039860 -0.522084 0.841280 
-0.198257 0.894806 -0.304826 
0.473945 0.834475 -0.252208 
0.632626 0.016126 0.705772 
0.030462 -0.199457 0.975695 
-0.112537 -0.193119 0.972142 
-0.373893 -0.377061 0.844031 
-0.376723 -0.130280 0.905450 
-0.802208 0.544817 0.123196 
-0.869841 -0.223571 -0.241933 
0.160888 0.910483 -0.333052 
0.007468 0.802257 -0.570580 
0.104069 0.367146 -0.917226 



0.301741 0.507113 -0.798909 
-0.023848 -0.005394 -0.995609 
0.158629 -0.006953 -0.985479 
0.065661 0.192695 -0.976218 
0.028183 0.094604 -0.992919 
0.091162 0.094314 -0.989423 

5 5 6 7 9 4 1 

8 17 16 14 13 12 11 19 10 2 

9 20 10 17 18 23 26 25 21 3 

6 20 19 11 2 22 4 
5 3 2 11 12 1 5 

5 9 4 13 12 1 6 

4 5 14 13 4 7 

5 6 15 16 14 5 8 

4 7 8 15 6 9 

5 15 16 17 18 8 10 

S 3 24 23 18 8 7 9 1 11 

7 3 24 27 25 21 22 2 12 
3 19 20 10 13 

3 21 22 14 

4 24 27 26 23 15 
3 26 27 25 16 



Something is lacking when viewing the polyhedra displayed in the two-dimensional images shown in 
this paper. To alleviate this difficulty, the Matlab and OFF files associated to each polyhedron that 
is constructed in this paper are included as supplementary materials on the website of Experimental 
Mathematics. See the website for full details on the use of Geomview. 



2.2 The desired polyhedron as a solution to 4A^ quadratic equations in 47V 
unknowns 

The proof of Andreev's Theorem gives that ac '■ Vq' Ac is a homeomorphism, so the problem of 
constructing a polyhedron P realizing (C, a) can be expressed as the problem of finding a solution P to 
the equation ac{P) = a. 

Instead of working in P^, we write the desired polyhedron as a solution of a system of 4A^ quadratic 
equations in AN variables, where is the number of faces of C. Our solution is vectors ■ ■ ■ , vm £ 
E^'^ satisfying 



• (Vi.Vi) = 1 and 



• {vj, Vj) ~ ~ cos(aij) if faces i and j are adjacent in C and their common edge is assigned dihedral 
angle a^j . 

These equations impose E + N conditions on AN variables, where C has N faces and E edges. 

As mentioned in Section 11.31 we have E = 3N — 6, so we have imposed AN — 6 conditions on 
AN variables. We impose 6 additional conditions in order to have the same number of equations and 
unknowns. We normalize by requiring that a chosen vector Vi perpendicular to one of the faces agree 
with some given v (where v is chosen so that {v,v) = 1.) We then require that one of the vertices on the 
face perpendicular to is at a given point lu in the plane Py and that a vertex adjacent to this vertex be 
on a given line I in P-^ through w. One can check that these normalizations provide 3,2, and 1 additional 
equations respectively. (Notice that the six equations for this normalization are each linear.) We denote 
the normalization by a triple {v, w, I). 

We denote the resulting quadratic map by ^c,(ii,to,0 ■ I^^^ ~^ . Typically we will only mention 
the normalization when necessary. We denote the conditions described above for the right hand side of 
the equations F(x) ~ y a,s (a, 0), where the a from this pair is shorthand for the conditions {vi,Vi) = 1 
and {vj,Vj) = — cos(aij) if faces i and j are adjacent in C, and the represents the fact that the 
normalization (d, w, I) is satisfied. 

Andreev's Theorem asserts that if a G Ac, there is a real solution to Fc.(v,w,i)ix) = (a, 0) correspond- 
ing to N vectors vi, • • • , vat in E^'^ so that P = nr=o realizes the pair (C,a). 

There are many sensible ways to numerically solve a system of quadratic equations in the same 
number of equations as unknowns. These include the pre-packaged non-linear solvers in Matlab, Maple, 
and Mathematica, Newton's Method, as well as Groebner basis techniques and fancier quadratically 
constrained solvers. 

The difficulty is that with AN quadratic equations in AN unknowns, Bezout's Theorem states that 
there will typically be 16N^ solutions. On their own, these solvers cannot easily be adapted to find the 
specific solution corresponding to a convex polyhedron, without first finding all solutions (or at least all 
real solutions) and then examining each solution to check if it corresponds to the desired polyhedron. 
Since some solutions may be much harder to find than others, one could spend significant computation 
times pursuing solutions that aren't of interest. 

One way to ensure that the solution does correspond to a compact convex polyhedron is to use an 
iterative method, like Newton's Method, for which an initial condition that is sufficiently close to a given 
solution is garunteed to converge to that root, in combination with a homotopy that garuntees that the 
nearest root is always the root that corresponds to a compact convex polyhedron. This is our approach, 
which we describe in greater detail in the next few sections of the paper. We are not entirely sure that 
this method is faster than finding all of the roots with "brute force" and then checking each solution to 
see if it is the desired one, but our approach has the additional benefit that it explicitly follows Andreev's 
proof of existence, providing insight into how this proof works for specific examples. 

2.3 Newton's Method and Homotopy methods 

Given two vector spaces V and W of the same dimension and a mapping F : V W, the associated 
Newton map Np : V V is given by the formula 

7VF(x)-x-[i?F(x)]-i(F(x)). (1) 

If the roots of F are non-degenerate, i.e. DF{ri) is invertible for each root n of F, then the roots of F 
corresponds bijectively to super-attracting fixed points of Np. 

Kantorovich's Theorem |27] gives a precise lower bound on the size of the basin of attraction for a 
root. 

Theorem 2.1 (Kantorovich's Theorem). Let ao be a point in R", U an open neighborhood of ao in 
M", and F : U ^ M" a differentiable mapping with [DF{ao)] invertible. 
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Let Uq be the open ball of radius |[D_F(ao)]~'^_F(ao)| centered at ai — NF{ai). If Uq C U and 
[DF{x)] satisfies the Lipshitz condition \\ DF{ui) — DF{u2) \\< ^^|ui — U2I for all Ui,U2 G Uq, and if 
the inequality 

|F(ao)|-|[i^F(ao)]-ipM<i (2) 

is satisfied, then the equation F(x.) ~ has a unique solution in Uq, and Newton's Method with initial 
guess ao converges to it. 

For a proof of Kantorovich's Theorem see [23 or the original source |27) . 

While the dynamics near a fixed point can be easily understood by Kantorovich's Theorem, the global 
dynamics of Newton's Method can be very complicated, with loci of indeterminacy, and critical curves 
where DN is not injective. In fact, the dynamics of Newton's method to solve for the common roots of 
a pair of quadratic polynomials in is a field of active research j^SJ 201 • We expect that the global 
dynamics of the Newton map to solve fc,(u,tu,;)(-^) ~ (^'0) even significantly more complicated than 
those in |28l I4()|. In particular, we have no reason to expect that a general initial condition in i?^^ 
will converge under iteration of Np to any solution of Fc^(v,w,i)ix) ~ (a, 0) nor to the specific solution 
representing a convex compact polyhedron P. 

An approach that can sometimes be used to avoid the difficulties with the global dynamics of Newton's 
Method is the homotopy method. Suppose that you want to solve g{x) = y. The idea is to replace this 
equation by a family that depends continuously on a single variable: 

gt{xt) = yt 

so that gi is the same function as g and yi = y, while goix) = yo is an equation for which you already 
know a solution xq. 

Choose k points = ti,t2, ■ ■ - tk = 1. If fc is sufficiently large, then xt-^ uiay be in the basin of 
attraction of the Newton's method for /^^ (x) = y^^ . In this case, you can solve for Xt.^ and can attempt to 
solve for Xt^ using Newton's Method for ft^ (x) — yt^ with initial condition Xt2 ■ Repeating this procedure, 
if possible, leads to the solution xi — Xt^ . 

While this is obviously a very powerful method, there are many difficulties choosing appropriate paths 
gt{xt) = yt and appropriate subdivisions = ti, i2, ■ • ■ , ife = 1- It is necessary to check that the conditions 
for Kantorovich's Theorem are satisfied by xt^ for the equation ftj+i{x) — ytj+i- The biggest difficulty 
is to avoid the situation where the derivative ^/t is singular for some t. Such points are described as 
being in the discriminant variety and choosing paths that avoid the discriminant variety is a big program 
of research. These difficulties are discussed extensively by many authors including Shub and Smale in 

[inEiEziEniEHiiiniizi. 

The proof of Andreev's Theorem in |25[l41j provides an explicit path that we can use for a homotopy 
method to construct any simple polyhedron P as a continuous deformation of either the prism Pr^ or 
the split prism Dn, both of which can be easily constructed "by hand". We will use this path for our 
homotopy method: repeatedly using a polyhedron realizing a point on the path as initial condition and 
solving for a polyhedron slightly further on the path, chosen so that the dynamics of Newton's method 
converges to the correct solution of F. 

With a similar path we can use the homotopy method again to construct any truncated polyhedron 
for which Aq ^ 0. We take a continuous deformation of a simple polyhedron until the vertices to 
be truncated pass SooEI^, and then add a finite number of additional triangular faces intersecting the 
appropriate triples of faces perpendicularly. Compound polyhcdra are then constructed as gluings of a 
finite number of truncated polyhcdra. 

Proposition 2.2 The quadratic equation F has a uniform Lipshitz constant on R^^ depending only on 
the combinatorics C. 
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Proof: The proof is merely the observation that F is quadratic so each of the second derivatives are 
constant. □ 

While we have checked that F is Lipshitz, we make no effort to bound the norm of the derivative [DF] 
away from zero (hence avoiding the discriminant variety). In fact, for a typical problem this is very hard 
to do. Instead, we merely try the homotopy method with the path mentioned in the preceding paragraph 
and we show that the method works for all of the constructions that we attempt. It may be interesting 
to provide a more rigorous basis for our use of Newton's Method and the current choice of path. 

2.4 Deforming a given polyhedron using Newton's Method 

Given a polyhedron P realizing C with dihedral angles a G Aq, it is easy to use Newton's method to 
deform P into a new polyhedron P' having any other angles a' £ Ac- Since Ac is a convex polytope, 
choose the line segment between a and a' and subdivide this segment into K equally distributed points 
a = a°, a^, a^, • • ■ a"^^^^ = a'. Then we use Newton's method with initial condition corresponding to P 
to solve for a polyhedron Pi with dihedral angles a-'^. We then repeats, using Pi as initial condition 
for Newton's method to solve for a polyhedron P2 with dihedral angles a^, and continue until reaching 
P' realizing a'. If the homotopy method has worked, then each step of Newton's method converges; 
otherwise we can try a larger number of subdivisions A', or attempt to check if the path has hit the 
discriminant variety. 

In all of the calculations within this paper, when deforming the angles of a given polyhedron P within 
Ac, we use K between 100 and 300 subdivisions, although this is sometimes a significant overkill. 

We consider it sufficient to show how to use Newton's method to construct some polyhedron P with 
non-obtuse dihedral angles for every C that has Ac 7^ 0- From this P one can construct any other 
P' £ Vc using the deformation described above. (This ease with which one can deform the angles of a 
given polyhedron is an additional benefit of our homotopy method.) 

In the next sections we will see how to connect individual paths in Ac\ , • • ■ , Ac^. so as to construct 
compact polyhedra realizing Ci a a sequence of deformations of a compact polyhedron realizing Ck- 

2.5 Simple polyhedra and Whitehead moves 

Recall that if C is simple then (^, ■ • ■ , ^) G Ac- The goal of this section and the following is to 
demonstrate the construction of a polyhedron P realizing any simple C with these dihedral angles. 

Andreev's Theorem provides a sequence of elementary changes (Whitehead moves) to the reducing 
the combinatorics C to one of two the combinatorial polyhedra 01 Ptn depicted below. 



I 




I 



Prism with 10 faces Splitprism with 11 faces 

In this section we show how to create polyhedra realizing Djsi and Prj^ and how to do the Whitehead 
moves using Newton's method. 
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Lemma 2.3 Let Pr^ and Dn be the abstract polyhedra corresponding to the N-faced prism and the N- 
faced "split prism", as illustrated below. If N > A, Vp^,^ is nonempty and if N > 7, is nonempty. 

Construction: Construct a regular polygon with N — 2 sides in the disc model for H^. (TV — 2 > 3, 
since > 5.) We can do this with the angles arbitrarily small. Now view as the equatorial plane of 
H^, and consider the hyperbolic planes perpendicular to the equatorial plane containing the sides of the 
polygon. In Euclidean geometry these arc hemispheres with centers on the boundary of the equatorial 
disc. The dihedral angles of these planes arc the angles of the polygon. 

Consider two hyperbolic planes close to the equatorial plane, one slightly above and one slightly 
beneath, both perpendicular to the z-axis. These will intersect the previous planes at angles slightly 
smaller than 7r/2. The region defined by these N planes makes a hyperbolic polyhedron realizing the cell 
structure of the prism. Note that our construction completes the proof of ProDOsition ll.61 for the special 
case C = Prif, N >5. 




Figure 2: 

For N > 7 we will construct hy cutting it into two prisms each with iV — 1 faces and the dihedral 
angles shown below. 



7r/2 




tt/4 



Figure 3: 

These angles satisfy Andreev's conditions (1) - (5) so we can use Newton's method to deform the 
prism constructed in the previous paragraph to have these angles. Gluing this prism is to its mirror 
image, the edges labeled 7r/2 on the outside disappear as edges, and the edges labeled on the outside by 
7r/4 glue together becoming an edge with dihedral angle tt/2. Hence, we have constructed a polyhedron 
realizing Dn, assuming N > 7. Notice that when N < 7 the combinatorics of coincides with that of 
PrN. □ 

Assume that the two vertices incident at an edge e are trivalent. A Whitehead move Wh{e) on edge e is 
given by the local change of the abstract polyhedron described in the following diagram. The Whitehead 
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move in the dual complex is dashed. Often we will find it convenient to describe the Whitehead move 
entirely in terms of the dual complex, in which case we write Wh{f). 



The following lemma appears in |25j : 

Lemma 2.4 Let the abstract polyhedron C be obtained from the simple abstract polyhedron C by a 
Whitehead move Wh(e). Then if is non-empty so is V^,. 

The proof constructs a sequence of polyhedra realizing C with dihedral angles chosen so that the 
edge e converges to a single point at infinity. A carefully chosen small perturbation of this limiting 
configuration results in a compact polyhedron realizing C" with non-obtuse dihedral angles. 

Suppose that we have a polyhedron realizing C with all dihedral angles equal to and choose a 
small e > 0. To implement a Whitehead move using the computer, we assign dihedral angle e to the the 
edge e and dihedral angle ^ to the four edges sharing an endpoint with e. Leaving the dihedral angles 
of the remaining edges the same, the resulting set of angles is in Ac and hence we can use Newton's 
Method to deform P into a polyhedron Pi realizing C with these new angles. 

If e was chosen small enough, Pi will be in the basin of attraction for a polyhedron realizing C" with 
the edge e' replacing e, the dihedral angle at e' equal to e, and all other dihedral angles as in Pi. We 
call the resulting polyhedron P2. Since C" is simple we can deform P2 to have all dihedral angles 



The following diagram shows these four steps when doing a Whitehead move on one of the edges of 
the dodecahedron. Here and elsewhere in this paper we use e approximately (A smaller e may be 
necessary when constructing polyhedra with a very large number of faces.) 

If we can find a sequence of combinatorial Whitehead moves reducing a given simple abstract poly- 
hedron C to either Pr^ or via a sequence of simple abstract polyhedra C'l, - ■ ■ ,Cn, we can use 
Newton's method to perform this sequence of Whitehead moves in the reverse order, constructing geo- 
metric polyhedra that realize Cjv, Cn~i, • • • , Ci, and finally C. Before explaining why such a sequence 
always exists, we demonstrate this process for the dodecahedron. 

We can then use Newton's Method to do the reverse sequence of Whitehead moves Wh{8,ll), 
Wh{A, 11), Wh{l, 2), Wh{9, 11), Wh{2, 4), Wh{l, 6), Wh{7, 11), Wh{6, 9), Wh{l, 5), and Wh{l, 4) geo- 
metrically, constructing the dodecahedron from -D12. The following diagram shows this process. 

2.6 A Lemma on Whitehead Moves 

The following lemma from |Hj and |25j is necessary to prove Andreev's Theorem and for our construction 
of simple polyhedra doing Whitehead moves geometrically with Newton's Method. 

Lemma 2.5 Whitehead Sequence 

Let C be a simple abstract polyhedron on which is not a prism. If C has N > 7 faces, C can 
be simplified to Dn by a finite sequence of Whitehead moves such that all of the intermediate abstract 
polyhedra are simple. 




Whitehead move on edge e 



Figure 4: 



obtaining P'. 
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Figure 5: Whitehead move Wh{e) on edge e of the dodecahedron 
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Theorem 6 in Andreev's original paper contains our Lcmma l2.5l Andreev's original proof of Theorem 
6 provides an algorithm to produce the Whitehead moves needed for this lemma, but the algorithm 
contains a glitch . An error was detected when the author tried to implement it as a computer program. 
The algorithm failed on the first test, when given a dodecahedron. 

In the instead using Wh{6, 9) for the fifth Whitehead move of the sequence described in the previous 
section, Andreev's algorithm uses either Wh{2,6) or Wh{2,5). In both cases it produces an abstract 
polyhedron which had a prismatic 3-circuit, see below: 




Figure 8: 

In combination with the computer-implemented Whitehead move described in the previous section, 
the sequence of Whitehead moves given in the proof of Lemma gives us the path that we will use for 
our homotopy method when constructing simple polyhedra. We provide an outline of the proof here that 
is sufficient to describe the sequence of Whitehead moves. Those who wish to see a complete proof may 
refer to ESUH]. 

Outline of the Proof of Lemma 12.51 

We assume that C ^ Pr^ is a simple abstract polyhedron with N > 7 faces. We will construct a 
sequence of Whitehead moves that change C to D^, so that no intermediate complex has a prismatic 
3-circuit. 

Find a vertex Voc of C* which is connected to the greatest number of other vertices. We will call 
the link of Uoo, a cycle of k vertices and k edges, the outer-polygon. Most of the work is to show that 
we can do Whitehead moves to increase fc to A'^ — 3 without introducing any prismatic 3-circuits during 
the process. Once this is completed, it is easy to change the resulting complex to by additional 
Whitehead moves. 

Let us set up some notation. Draw the dual complex C* in the plane with the vertex at infinity, 
and the outer polygon P surrounding the remaining vertices and triangles. We call the vertices inside of 
P interior vertices. All of the edges inside of P which do not have an endpoint on P are called interior 
edges. 

Note that the graph of interior vertices and edges is connected, since C* is simple. An interior vertex 
which is connected to only one other interior vertex will be called an endpoint. 




. endpoint 



Throughout this proof we will draw P and the interior edges and vertices of C* in black. The 
connections between P and the interior vertices will be grey. Connections between P and Voo will be 
black, if shown at all. 
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The link of an interior vertex v intersects P in a number of components ■ ■ , F" (possibly n ~ 0.) 
See the above figure. We say that v is connected to P in these components. Notice that since C* is 
simple, an endpoint is always connected to P in exactly one such component. 

Move 1 Suppose that there is an interior vertex A of C* which is connected to P in exactly one com- 
ponent consisting of exactly two consecutive vertices Q and R. The Whitehead move Wh[QR) on C* 
increases the length of the outer polygon by one, and introduces no prismatic 3-circuit. 




Move 2 Suppose that there is an interior vertex A that is connected to P in a component consisting of 
M consecutive vertices Qi, - ■ ■ ,Qm of P ( and possibly other components) . 

(a) If A is not an endpoint and M > 2, the sequence of Whitehead moves Wh{AQM), ■ ■ ■ ,Wh{AQ3) 
results in a complex in which A is connected to the same component of P in only Qi and Q2- These 
moves leave P unchanged, and introduce no prismatic 3-circuit. 

Q1Q2Q3 Qm-1 Qm Q1Q2Q3 Qm-iQm 




DAE DA 



(b) If A is an endpoint and AI > 3, the sequence of Whitehead moves 

Wh{AQM)-, ■ ■ ■ 1 Wh{AQ4) results in a complex in which A is connected to the same component of P in 
only Qi,Q2, and Q3. These moves leave P unchanged and introduce no prismatic 3-circuits. 




Note: In both parts (1) and (2), each of the Whitehead moves Wh{AQM) transfers the connection 
between A and Qm to a connection between the neighboring interior vertex E and Qm-i. This is helpful 
in case 2 later. 
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Move 3 Suppose that there is an interior vertex A whose link contains two distinct vertices X and Y 
of P. Then there are Whitehead moves which eliminate any component in which A is connected to P, if 
that component does not contain X orY . P is unchanged, and no prismatic 3-circuits will be introduced. 




Here A is connected to P in four components containing six vertices. We can eliminate connections 
of A to all of the components except for the single-point components X and Y. 

The proof that this move does not introduce any new prismatic 3-circuit is rather technical and 
depends essentially on the fact that A is connected to P in at least two other vertices X and Y. Andreev 
describes a nearly identical process to move |31 in his paper |S] on pages 333-334. However, he merely 
assumes that A is connected to P in at least one component in addition to the components being 
eliminated. He does not require that A is connected to P in at least two vertices outside of the components 
being eliminated. Andreev then asserts: "It is readily seen that all of the polyhedra obtained in this way 
are simple..." In fact, the Whitehead move demonstrated in Figure creates a prismatic 3-circuit. 

Having assumed this stronger (and incorrect) version of move 3, the remainder of Andreev's proof is 
relatively easy. Unfortunately, the situation pictured in Figure 12.61 is not uncommon (as we will see in 
Case 3 below!) Restricted to the weaker hypotheses of Move|31we will have to work a little bit harder. 

Using Move 1, Move 2, and Move 3 we check that if the length of P is less than — 3, then there is 
a sequence of Whitehead moves that increases the length of P by one without introducing any prismatic 
3-circuits. 

Case 1: An interior vertex that is not an endpoint connects to P in a component with two or more 
vertices, and possibly in other components. 

Apply Move 13 decreasing this component to two vertices. We can then apply Move 13 eliminating 
any other components since this component contains two vertices. Finally, apply Move to increase the 
length of the outer polygon by 1. 

Case 2: An interior vertex that is an endpoint is connected to more than three vertices of P. 

We assume that each of the interior vertices that are not endpoints arc connected to P in components 
consisting of single vertices, otherwise we are in Case 1. 

Let A be the endpoint which is connected to more than three vertices of P. By Move El part (2), 
there is a Whitehead move that transfers one of these connections to the interior vertex E that is next 
to A. Now, one of the components in which E is connected to P has exactly two vertices. The vertex 
E is not an endpoint since fc < — 3 implies that there are at least three interior vertices. Once this is 
done, we can apply Case 1. 

Case 3: Each interior vertex that is an endpoint is connected to exactly three vertices of P and 
each interior vertex which is not an endpoint is connected to P in components each consisting of a single 
vertex. 

First, notice that if the interior vertices and edges form a line, the restriction on how interior vertices 
are connected to P results in the prism, contrary to the assumption that C is not the prism. However, 
there are many complexes satisfying the hypotheses of this case which have interior vertices and edges 
forming a graph more complicated than a line: 
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For such complexes we need a very special sequence of Whitehead moves to increase the length of P. 

Pick an interior vertex which is an endpoint and label it Ii. Denote by Pi, P2, and P3 the three vertices 
of P to which Ii connects. Ii will be connected to a sequence of interior vertices /2, /s, • ■ ■ /m, ™ ^ 2, 
with Im the first interior vertex in the sequence that is connected to more than two other interior 
vertices. Vertex must exist by the assumption that the interior vertices don't form a line segment, the 
configuration that we ruled out above. By hypothesis, I2, ■ ■ ■ , Im can only connect to P in components 
which each consist of a vertex, hence each must be connected to Pi and to P3. Similarly, there is an 
interior vertex (call it X) which connects both to and to Pi and another vertex Y which connects to 
Im and P3. Vertex /,„ may connect to other vertices of P and other interior vertices, as shown on the 
left side of the following diagram. 




Now we describe a sequence of Whitehead moves that can be used to connect Im to P in only Pi and 
P2. This will allow us to use Move H to increase the length of P by one. 

First, using Move 13 we can eliminate all possible connections of /,„ to P in places other than Pi and 
P3. Next, we do the move Wh{ImP3) so that connects to P only in Pi. 




Next, we must do the moves W/i(/,„_iPi),...,W/i(/iPi), in that order (see the figure below.) 
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After this sequence of Whitehead moves we obtain the first diagram below, with /,„ connected to P 
exactly at Pi and P2, so that we can apply Move ^ to increase the length of P by the move Wh{PiP2), 
below. 





This concludes Case 3. 

Since C* must belong to one of these cases, we have seen that if the length of P is less than — 3, 
we can do Whitehead moves to increase it to — 3 without creating prismatic 3-circuits. Hence we can 
reduce to the case of two interior vertices, both of which must be endpoints. Then we can apply Move 
12 part (b) to decrease the number of connections between one of these two interior vertices and P to 
exactly 3. The result is the complex D^, as shown to the right below. 
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□ 



2.7 Construction of a "difficult" simple polyhedron 

We illustrate the algorithm described in the previous section by constructing a hyperbolic polyhedron 
realizing an abstract polyhedron for which Case 3 from the proof of Lemma 12.51 is necessary, and hence 
is particularly difhcult. 
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(4) (5) (6) 




(7) (8) (9) 




(10) (11) (12) 

Case 3 from the proof of Lemma 12. 51 is done in sub-figures (l)-(7). 
Case 1 follows in sub-figures (7)-(9) and again (for a different edge of P) 
in sub-figures (9)-(12). 



23 




Continuing from the last page, sub-figures (22)-(29) are another 
instance of Case 1. The figure has been straightened out between 
sub- figures (29) and (30). A sequence of final Whitehead moves 

is done in sub-figures (30)-(33) so that one of the two interior 
vertices is only connected to three points on the outer polygon 
This reduces the complex to Dig. 



Following the Whitehead moves backwards from £'*g back to Rig we find the following sequence of 
Whitehead moves. 

Wh{6, 13), Wh{6,9), Wh(3,&), Wh{9,12), Wh{6, 15), Wh{6,16), Wh{6,17), Wh{6,7), Wh{6,8), iyh(6,4), Wh{8,18), 
Wh{7,9), Wh{9,U), Wh{9,15), Wh{3,18), Wh{7,8), Wh{4,18), Wh{3,8), Wh{8,9), Wh{5,18), Wh{i,9), Wh{6,9), 
Wh(2,18), Wh{5,6), Wh{l,18), Wh(l,13), Wh{l,3), iy/i(3,4), VKh{4,5), Wh{2,5). 

We did this sequence of Whitehead moves geometrically, using Newton's Method. The result, starting 
with Dig, and realizing Rig is shown in Figure I^TI Each polyhedron is displayed in the conformal ball 
model. 



2.8 Truncation of vertices 

We have seen an outline of how to construct simple polyhedra. We now show how to construct all 
truncated polyhedra, except for the triangular prism, which wc have already constructed in Section 12.51 

Lemma 2.6 If Ac ^ 0, then there are points in Ac arbitrarily close to 
(7r/3,7r/3,...,7r/3). 

Proof: Simply check that if a G Ac, then the entire straight line path to (7r/3,7r/3, . . . ,7r/3), excluding 
the final point is in Ac- □ 
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9k» 










Figure 9: Construction of Ris from Dis using 30 Whitehead moves. 
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Thus we can assume that a is arbitrarily close to (tt/S, tt/S, • ■ • , tt/S), because once we have a poly- 
hedron realizing C with non-obtuse dihedral angles, we can deform it to have any dihedral angles in Ac-, 
as described in Section VIM Specifically, choose some < i5 < and assume that each component of a 
is within 5 oi ^. 

Let C be the C with each of the triangular faces /*' replaced by a single vertex vf^ . (Or, if C is the 
truncated triangular prism, let C be the prism.) Let a be the angles from a corresponding to the edges 
from C that are in C and let (3 = (ai + 2(5, 3.2 + 2(5, ■ • •) . (If C is the prism, re-number the edges so that the 
three edges forming the prismatic cycle are the first three, and choose (i = (ai, a2, ai, £15 -1-2(5, + ■ ■ ■).) 

Note that 5 was chosen so that /3 € A^. Then, the straight line path a(i) joining /3 to a (parameterized 
by t G (0, 1)) will remain in A^ except that the sum of the dihedral angles of edges meeting at each of 
the vertices vf will decrease past tt at some time ti S (0, 1). 

In |25II4H the authors use the path a(t) to construct a sequence of polyhedra P = Pq, -Pi, ■ ■ ■ , -P/v-i = 
P where P realizes C, and Pi is obtained from Pi^i truncating the vertices that become ideal when t = ti. 
(The truncation is done using a small perturbation.) Realizing P proves that 7^ 0, as needed for the 
proof of Andreev's Theorem. 

Because the proof from |25l I41| gives us a priori knowledge that compact polyhedra exist realizing 
each of the intermediate combinatorial structures exist, we can use Newton's method to deform the planes 
forming P to realize the angles in the entire path a{t) without truncating each vertex once it meets doo^^- 
We can then solve independently for the planes corresponding to the triangles in C so that each intersects 
the three appropriate planes at the three appropriate angles. 

We illustrate this construction below for the following pair (C, a), which has three truncations labeled 

fr, /|^ and fr. 




The path a.{t) described above works with any truncated C. For many C, such as the current one, a 
much easier path a(t) can be found satisfying conditions (1,3,4) and (5) from Andreev's Theorem for C, 
but for which the sum of the dihedral angles of the edges meeting at each v^" decreases past tt at some 
ti G (0, 1). Such a path is sufficient for our construction. For the current construction, C is shown below 
with edges labeled according to an appropriate path a(<). (Notation: a(t) ~ j3 for the edges labeled with 
a single angle (3, whereas a(i) = 77(1 — t) + -ft for the edges [77, 7].) 
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We constructed a polyhedron P realizing the pair (C,a(0)) and used Newton's method to deform 
the faces so that the dihedral angles follow the path a(t). After obtaining a non-compact polyhedron Pi 
realizing angles a(l) we truncated the vertices , , v\^' . The final result is the polyhedron: 




which will be used in the next section to form part of a compound polyhedron. 
2.9 Constructing compound polyhedra 

Any compound polyhedron can be constructed by gluing together a finite number of truncated polyhedra. 
We illustrate this construction for the polyhedron shown in Section [T!^ 

In general, one cuts along every prismatic 3-circuit which does not correspond to a triangular face. 
Here there is one such circuit which is labeled 7. We cut along 7 obtaining two combinatorial polyhedra 
for which every prismatic 3-circuit corresponds to a triangular face: 
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Figure 10: 




In this case, the diagram on the left is that for the polyhedron that we constructed in the previous 
section. The diagram on the right is that of the truncated triangular prism, which can also be easily 
constructed. 

We require that the new triangular faces F and F obtained by cutting along 7 be perpendicular to 
each of the other faces that they intersect. Then each face angle equals the dihedral angle outside of F, 
or F, that leads to that vertex. Because we obtained the two diagrams by cutting the original diagram 
along 7^the dihedral angles on the edges leading to F and F are the same and we naturally obtain that 
F and F have the same face angles, but are mirror images of each other. 



29 




F F 

These two polyhcdra glue perfectly together to form a polyhedron realizing (C, a) as shown in the 
following figure. 




3 Applications to discrete groups and polyhedral orbifolds 

Let P be a finite volume hyperbolic polyhedron having dihedral angles each of which is a proper integer 
sub-multiple of n. It is a well-known application of the Poincarc Polyhedron Theorem |2L)| that the group 
generated by reflections in the faces P forms a discrete subgroup Fp of Isom{M'^). Such groups have 
been extensively studied, see |55j . and the references therein. 

Given such a discrete reflection group Tp, we denote the corresponding orbifold by Op = H'^/Fp. 
We will use the term polyhedral orbifolds to describe orbifolds obtained in this way. (Note: often in the 
literature, the term "polyhedral orbifold" is used to describe the oriented double cover H'^/Fp, where Fp 
is index two subgroup consisting of orientation preserving elements of Fp.) See Thurston |51l Chapter 
13] and Reni [221 for more details on polyhedral orbifolds. 

We use the computer program described in this paper to construct examples from three classes of 
polyhedral orbifolds, the Lambert cubes [2H1, the Lobell orbifolds |29l I52L 1^ . and an mysterious orb- 
ifold described by Mednykh and Vcsnin whose 16-fold cover is a"hyperelliptic" compact hyperbolic 
manifold. 

We output the generators each reflection group as elements of 5*0(3, 1) into SnapPea [SS], computing 
volumes and length spectra of these orbifolds. For details on how SnapPea calculates the length spectrum 
refer to 
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3.1 Construction of Lambert Cubes 

A Lambert cube is a compact polyhedron realizing the combinatorial type of a cube, with three non- 
coplanar edges chosen and assigned dihedral angles a, /9, and 7, and the remaining edges assigned dihedral 
angles ^. It is easy to verify that if < a, /3, 7 < ^ then, such an assignment of dihedral angles satisfies the 
hypotheses of Andreev's Theorem. The resulting polyhedron is called the (a,/3,7) Lambert Cube, which 
we will denote by Pa/j.-y- Thus, there are discrete reflection groups generated in the faces of a Lambert 
Cube when a = ^, /3 = ^, and 7=7 for integers p,q,r > 2. We denote the corresponding orbifold 
OLambcrt(p, 9, J') In the following table, we present volumes and the lengths of the shortest geodesies for 
a sampling of Lambert Cubes for small p, q, and r: 



OLambcrt(3, 3, 3) 


Computed Volume: 0.324423 


Theoretical Volume: 0.3244234492 


Short Geodesies 


6 ml 1.087070 
6 ml 1.400257 

3 ml 1.601733 + i • 2.765750 
6 ml 1.864162 

4 ml 2.174140 


3 ml 1.087070 + i cot 27r/3 

6 ml 1.400257 + i-TT 

3 ml 1.790480 + ^-0.762413 

3 ml 2.138622 

6 ml 2.199243 + i - 2.436822 


OLambcrt(3,4, 5) 


Computed Volume: 0.479079 


Theoretical Volume: 0.4790790206 


Short Geodesies 


2 ml 0.622685 

1 ml 0.622685 + i ■ 2.513274 

1 ml 0.883748 + i • 1.570797 

1 ml 1.123387 

1 ml 1.245371 


1 ml 0.622685 + i - 1.256637 
3 ml 0.883748 
1 ml 0.883748 + i - 3.141592 
1 ml 1.123387 + i - 3.141593 


OLambcrt(4,4,4) 


Computed Volume: 0.554152 


Theoretical Volume: 0.5382759501 


Short Geodesies 


2 ml 0.175240 
1 ml 0.175240 + i • 1.108797 
1 ml 0.175240 + ^-0.739198 
1 ml 0.175240 + i • 1.847996 
1 ml 0.175240 + i • 2.587194 
1 ml 0.350479 


1 ml 0.175240 + ^-0.369599 
1 ml 0.175240 + ^-0.739198 
1 ml 0.175240 + i - 1.478396 
1 ml 0.175240 + i - 2.217595 
1 ml 0.175240 + i - 2.956793 


OLambcrt(5, 8, 12) 


Computed Volume: 0.768801 


Theoretical Volume: 0.7688005863 


Short Geodesies 


3 ml 0.407809 
1 ml 0.407809 + i • 1.047198 
1 ml 0.407809 + i ■ 2.094396 
1 ml 0.407809 + ^-3.141592 
1 ml 0.643110 + i- 0.785398 
1 ml 0.643110 + i - 2.356194 


1 ml 0.407809 + i - 0.523599 
1 ml 0.407809 + i - 1.570797 

1 ml 0.407809 + i - 2.617995 

2 ml 0.643110 

1 ml 0.643110 + i - 1.570796 
1 ml 0.643110 + ^-3.141593 



The format of the lists of Geodesic lengths presented in this and in the following tables, is the same 
as that presented by SnapPea. The first entry is the multiplicity of distinct geodesies having the same 
complex length. The second entry is either "ml" to indicate that the geodesic has the topological type of 
a mirrored interval or is empty, if the geodesic has the topological type of a circle. The third entry is the 
complex length. Nearly all of the short geodesies that we present in these tables are mirrored intervals 
because our orbifolds are mirrored polyhedra and because we have only listed rather short geodesies. 

Also notice, that while SnapPea provides many more digits of precision for the geodesic length, we 
have rounded to the first 6 decimal places in order to group geodesies that are likely to correspond to the 
same class, but weren't listed that way due to numerical imprecision. 

The volumes of Lambert cubes have been explicitly calculated by R. Kellerhals |2H1- If we write 
'^iViO ^ ^{v+O^-^iv^O^ where A is the well-known Lobachevskii function A(a;) = — Jq log |2 sin(t)|(it, 
then 
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Vo\{Pc.,0.,) = - (A(a, e) + A(/3, 6) + A(7, 0) - 2 ■ A (|, - A(0, 9)) 



(3) 



where 6, with < 6* < ^ is the parameter defined by: 

tan^ (61) = p+ ^p^ + L^APN^, 

+ AP + 1 

n = . and 

^ 2 

L = tana, M = tan/3, iV = tan7. 

The column in the above table labeled "approximate volume" gives the volume of OLambcrt(P: 97 '') 
computed using SnapPea, while the column labeled "actual volume" gives the volume of OLambcrt(p, <7, r) 
computed using Equation |3| 

3.2 Construction of Lobell Orbifolds 

For each n > 5, there is an radially symmetric combinatorial polyhedron having two n-sided faces 
and having 2n faces with 5 sides which provides a natural generalization of the dodecahedron. This 
combinatorial polyhedron is depicted below for n = 8. 



Andrcev's Theorem provides the existence of a compact right angled polyhedron i?„ realizing this 
abstract polyhedron because it contains no prismatic 3-circuits or prismatic 4-circuits. (In fact, the work 
of Lobell predates Andreev by many years, and one can also verify the existence of Rn as an appropriate 
truncation and gluing of tetrahcdra.) We refer to the group generated by reflections in the faces of i?„ 
by r„ and the corresponding orbifold OL6bou(?^) = IHI'^/r„. 

Historical Note: 

While we restrict our attention to the orbifold Oi^oheiiin) in this paper, the reader may wish to notice 
that the first example of a closed hyperbolic manifolds was constructed by Lobell [221 in 1931 by an 
appropriate gluing of 8 copies of R5. Generalizing this notion, Vesnin |52j has described a convenient 
algebraic method to construct a torsion free subgroup C r„ of index 8. This, in the n-th Lobell 
Manifold is the compact, orientable, hyperbolic manifold AfL6bcii("-) := H-^/Fj^. Naturally, ML6bcii(«) 
is an 8-fold (orbifold) cover of OL6bGii(?^)- We refer the reader to the nice exposition in jS^lEll for the 
details. The delightful paper by Reni |36| provides further details on the construction of hyperbolic 
manifolds and orbifolds and finite covers of right-angled polyhedra. 

We include the following table of data computed using SnapPea for the n = 5, • ■ • , 8 Lobell orbifolds. 
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OL6bcii(5) (Dodecahedron) 


Computed Vohime: 4.306208 


Theoretical Volume: 4.3062076007 


Short Geodesies 


60 ml 2.122550 

60 ml 2.938703 

126 ml 3.233843 

60 ml 3.783112 + i ■ 1.376928 

12 3.835986 + i ■ TT 

60 ml 3.966774 


60 ml 2.122550 + i • tt 

60 ml 2.938703 + i ■ w 

60 ml 3.579641 

12 3.835986 

60 ml 3.835986 + i • tt 

60 ml 4.0270318 + i ■ 2.264758 


OL6bcll(6) 


Computed Vohime: 6.023046 


Theoretical Volume: 6.0230460200 


Short Geodesies 


36 ml 1.762747 

37 ml 2.292431 
48 ml 2.633916 
36 ml 2.887271 
48 ml 3.088970 
24 ml 3.256614 


12 ml 1.762747 + i-TT 

12 ml 2.292431 + i • TT 

24 ml 2.633916 + i • tt 

24 ml 2.887271 + i • tt 

12 ml 3.154720 + i • 1.312496 

36 ml 3.256614 + i • tt 


OLobell(7) 


Computed Volume: 7.563249 


Theoretical Volume: 7.5632490914 


Short Geodesies 


42 ml 1.611051 

1 ml 1.823106 

14 ml 2.388409 + i ■ ir 

14 ml 2.512394 + i • tt 

70 ml 2.898149 

42 ml 2.898149 + i • tt 


14 ml 1.611051 + i-7r 

42 ml 2.388409 

14 ml 2.512394 

14 ml 2.601666 

14 ml 2.898149 + i ■ 1.280529 

14 mi 3.031090 + i • TT 


OL6bcll(8) 


Computed Volume: 9.019053 


Theoretical Volume: 9.0190527274 


Short Geodesies 


49 ml 1.528571 

80 ml 2.448452 

16 ml 2.760884 + i ■ 1.261789 

48 ml 2.914035 + i ■ tt 

32 ml 3.057142 + i ■ tt 

64 ml 3.553688 


16 ml 1.528571 + i-7r 

32 ml 2.448452 + i • tt 

32 ml 2.914035 

160 ml 3.057142 

16 ml 3.461816 + i ■ 2.650944 

32 ml 3.553688 + i ■ tt 



The column labeled "Computed Volume" gives the volume as computed in SnapPea, whereas "Theoretical 
Volume" provides the volume of OLobciil"-) using explicit formula from |53| . (In fact we have divided by 
8 the volume formula presented in [HSI, because they study the volume of the 8-fold cover A/L6bcn('T^)-) If 
we let = f - arccos [ ^cosl^/n) ) ' ^hen 

Vol(OLobcii(n)) = I (2A{9) + A + ^) + A - ^) - A (20 + I) ) . (4) 

where A is the Lobachevskii function. 

Notice that for each of the Lobell orbifolds that we computed, the volume computed in SnapPea 
agrees perfectly (within the six digits of precision available) with that given by Equation^ 

3.3 An orbifold due to Mednykh and Vesnin 

In a very similar way to the construction of Lobell manifolds, Mednykh and Vesnin describe in |35| a 
compact three-dimensional hyperbolic manifold G which forms a 2-fold branched cover over the sphere 
S'^. They call manifolds with such a covering property over S'^ "hyperelliptic," generalizing the classical 
notion of hyperelliptic Riemann surfaces. See also |32l 1341 

The combinatorial polyhedron considered by Mednykh and Vesnin (and apparently originally due to 
Grinbergs) is depicted below. 
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This abstract polyhedron has no prismatic 3-circuits or prismatic 4-circuits, so Andrccv's Theorem 
garuntees the existence of a polyhedron Rmv realizing it with tt/2 dihedral angles. We denote the group 
generated by reflections in the faces of Rmv by Tmv and the orbifold by Omv- Combinatorial details 
on the construction of Mmv as a 16-fold cover of Omv can be found in p?5] . 

The following table contains invariants of the orbifold Omv = ^^/^mv obtained by entering an 
explicit list of generators for Tmv into SnapPea. 



Oaiv 


Computed Volume: 6.023046 


Theoretical Volume: unknown 


Short Geodesies 


9 ml 0.989308 


3 ml 0.989308 + i ■ it 




9 ml 1.183451 


3 ml 1.183451 + i -TT 




18 ml 1.834468 


6 ml 1.834468 + i-TT 




18 ml 1.859890 


6 ml 1.859890 + i-TT 




27 ml 1.882318 


9 ml 1.882318 + i-TT 




6 ml 1.978616 


3 ml 1.978616 + i-TT 




9 ml 2.214787 


3 ml 2.214787 + i-TT 




18 ml 2.252719 


6 ml 2.252719 + i-TT 




6 ml 2.366902 


3 ml 2.366902 + i • TT 




6 ml 2.433170 


6 ml 2.433170 + i • TT 




6 ml 2.446977 


6 ml 2.446977 + i-TT 



As an application, we obtain the estimate VoI{Mmv) = 16 • 15.608119 = 249.729904 using that Mmv 
is a 16-fold orbifold cover over Omv- 

3.4 Spectral Staircases 

For a given hyperbolic manifold or orbifold M, the "spectral staircase" is a plot of the number of closed 
geodesies of length less than I, which we denote N{1), versus I. (In fact, it is much more common to plot 
log(A^(/)) due to the exponential grown predicted by ^ below.) The spectral staircase provides both a 
nice way to graphically display the spectrum of M and a illustration of the classical result of Margulis 
|31j . who proves the following universal formula for the asymptotics of N{1): 

N{1) ~ as / oo (5) 

Tt 

where the constant r is the topological entropy, which for hyperbolic space H'^ is given by r = d — 1. 
For an exposition and nice experimental works considering spectral staircases, see and the references 
within. 
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Figure 11: Spectral Staircases for the OLambort(3, 4, 5) (solid line), OL6bcii(6) (dotted line), and Omv 
(dashed line). 

We compute these spectral staircases for OLambcrt(3, 4, 5), OL6bcii(6), and Omv displaying the results 
in Figure 1X11 (The data for OLambcit(3, 4, 5) ends at roughly / — 3.8. SnapPea encounters an error 
computing at this length, probably due to the comparatively small dihedral angles of 0Lambcrt(3, 4, 5).) 

4 Questions for further study 

We present a non-comprehensive list of interesting questions for further study: 

1 Determine if there is a faster way of computing Andrecv Polyhcdra. 

lb Related question: determine if CirclePack j^O] can be used to construct compact hyperbolic poly- 
hedra and their reflection groups. If so, it would provide a faster method of construction. 

2 Construct manifold covers of the polyhedral orbifolds that were considered in Section O including 
the Lobell Manifolds the "Small Covers of the Dodecahedron" ^7], and the Hyperelliptic 
Manifold Such a construction, could potentially lead to computations of many additional 
interesting invariants of these manifolds using SnapPea, as well as on drilling and Dehn fillings on 
them (which would also be possible in SnapPea). 
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2b Related question: use the program SNAP ^S] to compute arithmetic invariants for these manifolds. 

2c Related question: using SNAP, or the ideas used in SNAP [TJ, study the arithmetic invariants of 
polyhedral reflection groups. 

3 Perform a study of volumes of hyperbolic polyhedra corresponding to general angles in Ac- (While 
SnapPea computes volumes only for polyhedra with discrete reflection groups, the functions from 
the SnapPea kernel could probably be used for this more general study.) 



Appendix: 

A Using the program to construct compact polyhedra and their 
reflection groups 

The computer program described in this paper is a functional but slightly rough collection of Matlab 
(or Octave scripts. Using a single command the program produces a polyhedron realizing a simple 
abstract polyhedron C with all dihedral angles However, one must do a little bit more work to 
construct a polyhedron realizing (C, a) of truncated or compound type. These later steps are not as 
automatic in the program, but one can follow the description in this paper step by step to do them "by 
hand." Please see the README file enclosed with the program for further information. 

There are two ways to output a polyhedron that has been constructed using this program: the Object 
File Format (e.g. "filename. off") and the Generators File Format (e.g. "filename.txt"). The Object File 
Format output can be read into Geomview [3] and displayed nicely in the Poincare Ball model there. The 
Generators File Format is for SnapPea |58) and (if the polyhedron output this way has dihedral angles 
that are proper integer sub-multiples of tt) this file can be loaded into SnapPea. 

An analysis of the computational complexity of this method would be quite involved and is not feasible 
at this point. In fact, the computational complexity of Newton's Method is quite difficult |47ll46[H51l49l l7|. 

In practice, on a contemporary PC running Linux, the most complicated construction in this paper 
took approximately 2 minutes; it is by no means fast, but certainly usable. We expect that the program 
described in this paper can be used to construct all polyhedra having no more than 30 faces and having 
dihedral angles bounded away from dAc (with the exception of the part of dAc corresponding to dihedral 
angles 7r/2, where there should be no problem.) Certainly the number of iteration steps in each homotopy 
(the parameter K from Section and the parameter e used in the Whitehead move (see Section 
may need to be modified in special circumstances. With improvements of the program, perhaps by 
implementing it in a faster programming language, and, if necessary, in a higher precision arithmetic, we 
expect the same to be possible for up to 100 faces, or possibly more. 
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